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Abstract 



If tan f5 is large, down-type quark mass matrices and Yukawa couplings cannot 
be simultaneously diagonalized, and flavour violating couplings of the neutral 
Higgs bosons are induced at the 1-loop level. These couplings lead to Higgs- 
mediated contributions to the decays Bs /"^/^~ and Bd t^t~, at a 
level that might be of interest for the current Tevatron run, or possibly, at 
-B-factories. We evaluate the branching ratios for these decays within the 
framework of minimal gravity-, gauge- and anomaly- mediated SUSY breaking 
models, and also in SU{5) supergravity models with non-universal gaugino 
mass parameters at the GUT scale. We find that the contribution from gluino 
loops, which seems to have been left out in recent phenomenological analyses, 
is significant. We explore how the branching fraction varies in these models, 
emphasizing parameter regions consistent with other observations. 
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I. INTRODUCTION 



Supersymmetry provides a promising way to stabilize the electroweak scale if superpart- 
ners are lighter than 0{1) TeV [1]. All SUSY theories [2] necessarily contain many scalar 
fields, resulting in multiple potential sources of fiavour violation. Indeed, in constructing 
phenomenologically viable models, care has to be taken to make sure that fiavour violation 
is sufficiently suppressed. Within the framework of the Minimal Siipersymmetric Standard 
Model (MSSM) with conserved ii-parity, this is ensured (at tree level) by requiring that the 
matter supermultiplets with weak isospin T3 = 1/2 couple only to the Higgs superfield h^, 
while those with — —1/2 couple just to the Higgs superfield ha. 

At the one loop level, however, a coupling of hu to down type fermions is induced. This 
induced coupling leads to a new contribution [3], proportional to (hu), to the down type 
fermion mass matrix. Although this contribution is suppressed by a loop factor relative to 
the tree-level contribution, this suppression is partially compensated if the ratio = tan /3 
is sufficiently large. As a result, down type Yukawa coupling matrices and down type quark 
mass matrices are no longer diagonalized by the same transformation, and flavour violating 
couplings of neutral Higgs scalars h, H and A emerge. Of course, in the limit of large rriA, 
the Higgs sector becomes equivalent to the Standard Model (SM) Higgs sector with a light 
Higgs boson h ~ hsM, and the effects of fiavour violation decouple from the low energy 
theory . The interesting feature is that the fiavour-violating couphngs of h, H and A, do 
not decouple for large superparticle mass parameters: being dimensionless, these couplings 
depend only on ratios of these mass parameters^, and so, remain finite even for very large 
values of SUSY mass parameters. 

As pointed out by many authors [4-14], this fiavour violating neutral Higgs boson cou- 
pling results in a potentially observable branching fraction for the decay Bs — > A*"*"/^" medi- 
ated by the neutral Higgs bosons, h, H and A. If tan is sufficiently large the most important 
contribution to the amplitude for this decay then scales as tan^ (3: two of these factors arise 
via the down type quark and lepton Yukawa couplings which are each oc 1/ cos/? ^ tan/?, 
while the last factor arises because the offset between the down quark Yukawa coupling 
matrices and the mass matrices discussed above also increases with tan/3. 

Within the SM, the branching fraction for this decay is ~ 3.4 x 10~^ [15]. This is very far 
from its current experimental upper limit of 2.6 x 10^^ obtained by the CDF collaboration 
[16]. Within the SUSY framework, however, this branching fraction can be enhanced by a 
very large factor and the SUSY contribution may dominate its SM counterpart if tan/S is 
sufficiently large. Indeed, the CDF bound on this branching fraction already constrains some 
portions of the SUSY parameter space [11,10,17]. It is, therefore, reasonable to expect that 
the data from the Fermilab Main Injector will probe significant ranges of SUSY parameters. 
Moreover, because the SUSY amplitude does not decouple for large sparticle masses, it is 
possible that the CDF or D0 experiments may detect new physics via a measurement of 
B{Bs — > A*'*'//") even without direct detection of sparticles at the Tevatron. 

Many of the analyses [6-9,12,13] of Bg — > have been performed within the MSSM 



These not only include sparticle masses, but also the superpotential parameter and also the 
soft SUSY breaking 74-parameters. 
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framework. Since sparticle mass matrices as well as soft SUSY breaking scalar coupling 
matrices both include intrinsic sources of flavour violation, strictly speaking the framework 
is too broad to be predictive.^ Even with assumptions that set tree-level flavour violating 
effects to zero, the model still has too many free parameters to correlate the branching 
fraction for Bg decay with other observables. 

This situation is in sharp contrast to that in highly predictive scenarios such as the 
minimal supergravity (mSUGRA) model [18], the gauge-mediated SUSY breaking (GMSB) 
model [19], or the minimal anomaly-mediated SUSY breaking (mAMSB) model [20] in which 
much of SUSY phenomenology is analysed these days. In each of these scenarios, SUSY is 
assumed to be dynamically broken in a "hidden sector", comprised of fields that interact 
with SM particles (and their superpartners) only via gravity. It is the mechanism by which 
SUSY breaking is communicated to the superpartners of SM particles that distinguishes 
these frameworks from one another. These models are all completely specified by just a 
handful of parameters, usually defined at a scale, much higher than the weak scale, where 
the physics is simple. 

The main goal of this paper is to examine the prospects for observing the decay 
Bs — > fi'^fi~ in Tevatron experiments within these framework of well-studied supersym- 
metry models. We also consider gravity-mediated SUSY breaking models based on SU{5) 
where the order parameter for SUSY breaking also breaks the GUT symmetry [21], leading 
to non-universal GUT scale gaugino mass parameters [22]. With an integrated luminosity 
of 2 fb~^, experiments at the Tevatron are expected to be sensitive to a branching fraction 
below ~ 10"^. With a still bigger data sample (that is expected to accummulate before the 
Large Hadron Collider (LHC) begins operation) the sensitivity should be even greater. 

Of course, while detection of a signal would point to new physics, it would not necessarily 
establish this new physics to be supersymmetry. Within these specific scenarios, however, it 
is possible to predict other signals that might also be simultaneously present if an observed 
Bs fi~^fi~ signal is to be attributed to any particular supersymmetric framework. The first 
examination of the Higgs-mediated Bg fi^fjr decay in the supersymmetric context that 
we know about was performed within the mSUGRA framework [4]. However, it is only in the 
past year that the observabihty of this signal as a function of mSUGRA parameters been 
systematically investigated [10,11]. A corresponding study [14] of this signal within the 
gauge-mediated and the anomaly-mediated SUSY breaking scenarios also appeared while 
this paper was in preparation. As far as we are aware, there is no study of this decay 
in models with non-universal gaugino masses. As discussed in the next section we have 
improved upon the existing phenomenological analyses, which consider fiavour violating 
effects just from chargino loops, in that we also include effects from gluino loops which we 
find to be substantial. 

The Higgs-mediated mechanism for the ;U-pair decay of Bs mesons also allows the decay 
Bg — > T^T~ . In fact, since the leptons couple via their Yukawa interactions, the branching 
ratio for the latter decays would be expected to be enhanced by (m^/m^)^. Unfortunately, 
identification of r pair decays of Bg would be very difficult in Tevatron experiments because 



^It would, of course, be possible to obtain constraints on flavour-mixing elements of the squark 
mass matrix or 74-parameters. 
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their invariant mass cannot be accurately reconstructed. Moreover, Bg pair production is 
not kinematically accessible at the i?-factories currently in operation at SLAC or KEK. 
The BELLE and BABAR experiments at these facilities have already accumulated ~ 70M 
Ba mesons, and this sample is expected to increase to ~ 500 — lOOOM B^s in a few years. 
These considerations motivated us to also examine the branching fraction for B^ — > t'^t~. 
Although this decay is suppressed by Kobayashi-Maskawa (KM) matrix elements relative to 
corresponding decays of Bs, this is offset by the larger r Yukawa coupling to the Higgs bosons. 
We are not aware of any studies of the sensitivity of BELLE or BABAR to B^ —>■ t~^t~ 
decays. We will, therefore, confine ourselves to mapping out the branching fraction for this 
decay without making any representation of its detectability in experiments at S-factories. 

The same flavor-violating couplings of neutral Higgs bosons that we have discussed can 
also lead to significant contributions to other processes, such as the exclusive decays Bs — > 
Ki~^i~ and Bg £'^i~j, the semi-inclusive decay Bg — > Xsi~^i~, and to AF = 2 processes 
like Bd g and X-meson mixing [23] that may be probed by experiment. Their leptonic cousins 
lead to T — > ii/iiJ, and — > eee decays with branching fractions that might potentially be 
within reach of future experiments [24] . 

The remainder of this paper is organized as follows. In Sec. II we describe our computa- 
tion of the amplitude for the decay, focussing on the improvements that we have made, and 
on the approximations that we have used to simplify the analysis. Sec. Ill contains our main 
results. Here we analyse the branching fractions for the decays Bg — > and B^ — > t'^t~ 

within the mSUGRA, mGMSB and mAMSB models, as well as in SU{5) supergravity mod- 
els with non-universal GUT scale gaugino masses. We also comment on other observables 
(with emphasis on the other flavour-violating decay b s'~f) for parameter regions where 
B{Bg — > may be observable at the Tevatron. We conclude in Sec. IV with a summary 

of our results. 

II. HIGGS-MEDIATED LEPTONIC DECAYS OF Bgd MESONS 

At tree level, the down (up) type quarks couple to the Higgs fleld (/i^) via the Yukawa 
coupling matrix (fu). Flavour violation in Higgs fleld interactions occurs because couplings 
of the fleld to down type quarks is induced at the one loop level. There are two distinct 
sets of SUSY sources of this coupling at the 1-loop level [5] . Down type quarks can couple 
to hi via 

1. a chargino up-type squark loop, with the Higgs fleld attaching to the squarks via the 
Afu elements of the soft SUSY breaking trilinear scalar couphng matrix, or 

2. via a gluino down- type squark loop, with the Higgs fleld attaching to the squark via 
supersymmetric interactions proportional to the matrix /xfd. We ignore analogous 
contributions from neutralino loops; see however, Ref. [12]. 

The gluino loop, at flrst sight, appears to give no flavour violation. In the flavour basis, 
gluino interactions are flavour diagonal, so that the Higgs-squark-squark coupling is the only 
source of fiavour violation. Since the fiavour structure of this coupling is proportional to 
the down quark Yukawa matrix (i.e. the same as the tree-level Yukawa matrix) it does 
not lead to fiavour violation. This reasoning is incorrect because the down-squark matrices 
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need not be flavour diagonal. In other words, the gluino contribution to flavour violation 
depends on the a priori unspecified structure of the down squark mass matrices. The flavour- 
violating coupling from the chargino contribution above depends on the essentially unknown 
matrix Afu- These considerations make it clear why it is difficult to make predictions of 
the flavour violating branching fractions without resort to specific models, although as we 
said it is possible to constrain certain (flavour-mixing) matrix elements in a more general 
analysis. We will, in the following, focus on the four models introduced in Sec. I. The reader 
should, however, keep in mind that our results for the branching fraction are specific to these 
models, and that simple modifications to the flavour structure in the sparticle sector could 
lead to a very different answer. 

It is instructive to analyse the flavour structure of the gluino- induced flavour violation. 
Following Ref. [5], we will work in the quark basis where the up type Yukawa coupling matrix 
is diagonal, and the down quark Yukawa coupling matrix has the form, = DV^j^, where 
D is the matrix obtained by diagonalizing and Vkm is the Kobayashi-Maskawa matrix 
(uncorrected for SUSY contributions). In the mass basis for quarks, the flavour 

structure of the gluino contribution can be written as, 

(ZR)ji(Z]jDZi)zfe(Zi)fciCo(m|^^, m|^^, m?), (1) 

where m| (m| ) is the kth (/th) eigenvalue of the squared mass matrix for left (right) 
type squarks, Z]^ {Zr) is a unitary matrix that transforms the superpartners of (d^) to 
the basis in which the left (right) down squark squared mass matrix is diagonal and Co is a 
loop function defined by, 

.2 /„,2x ..,2 /^2xn 



Co{x,y,z) = — 

— 



y 



r^2 _ y2 



\x^J — z^ ) 



(2) 



The following features are worth noting: 



If the different types of left and right down squarks have common masses rxiL and mR 
(not necessarily equal), the function Co will be independent of k and I and the flavour 
violating contribution from (1) vanishes. 

In a wide class of models the left (and, separately, the right) squark mass squared 
matrices are proportional to the unit matrix at some scale, and only Yukawa inter- 
actions distinguish between flavours.^ Then, the matrix is diagonal in the basis 
of superpartners of d^, and = I. This may be seen by noting that the 1-loop 
renormalization group equation for the singlet down squark mass matrix depends only 
on the down type Yukawa coupling matrix. The same is not the case for the doublet 
squarks since the renormalization group equation now depends also on the up type 
Yukawa couplings. In all such models, gluino-induced flavour violation depends only 
upon the mass splittings in the left squark sector. Note that all the models that we 
consider all fall in this category. 



^Operationally, this means that the squark mass matrices have the form, m? 



1 + CLfdfd + c'^fufu and m? 



dR 



m 



R 



1 -|- csfdf^ , with cl, c'j^ and cr as constants. 



dL 
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A. Flavour-violating Higgs boson interactions 



In our analysis we follow the approach in Ref. [5] and develop an effective Lagrangian 
to describe the induced flavour violation in the neutral Higgs boson sector. We will assume 
that squarks with the same quantum numbers have a common mass at some energy scale, 
and further that superpotential Yukawa interactions are the sole source of flavour violation 
in the remainder of this analysis. We take the SU (3) gaugino mass parameter to be positive 
by convention: other gaugino mass parameters may have either sign. As in Ref. [5], we work 
in the basis where fu is diagonal and neglect masses for d and s quarks. Neglecting terms 
that are second order in the Wolfenstein parameter A, the effective couplings of down quarks 
to the Higgs fields of the MSSM can be written as,^ 

-£e// = DRi^QLhd + DRij,[agM'^ + a^M'^iu + a^M'']QLhl + h.c. , (3) 

with 
and 

= diag(Co^j, C^,:l Cl%), Cf;] ^ Co{mg, m^^^^,m^J , 

= diag(0, 0, CI,), CI, ^ Co{^^, mi^,mij , (5) 

- diag(Co-i, Co's), Cl^j = Co{M2, m,^,_ J . 

The loop function Co{x,y,z) is given in Eq. (2). Note that the right squark index in the 
matrix Ai^ in (3) is fixed by the corresponding quark; matrix multiplication is implied with 
the left squark index of the "diagonal matrix" A4^. 

The second term of (3) contains the loop-induced couphng of /i° to down quarks. The 
first entry in this term comes from the gluino loop, while the last two entries arise due to 
chargino loops. We see immediately that if Qj^l have a common mass, there is no flavour 
violation from the gluino term since it then has exactly the same flavour structure as the tree- 
level term. The term depending on A^" is proportional to the squared top Yukawa coupling, 
and arises from the Higgsino components of the chargino, with the fermion chirality flip 
proportional to the Higgsino mass fi. The last term in (3), which comes from gaugino- 
higgsino mixing and again vanishes if left type squarks have a common mass. It is, generally 
speaking, much smaller than the other two terms. 

The flavour changing coupling between second and third generation down quarks is given 

by, 

fb 

C^FCNC = -^—rV^^ytsXFchRSL{cosl3hl - sin phd) 
sm p 

fb 

ytb^tsXFcbRSL[cos{P + a)h - sin(/? + a)H - lA] , (6) 



V2sin/3 



^This is the analogue of Eq. (7) of Ref. [5]. 
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where a is the mixing angle in the neutral scalar Higgs boson sector, is the "physical h 
Yukawa coupling" defined by [5] , 

/, = + (a,Co^,| + + a^Cl^f^) tan /3] , (7) 



and 



[%(C'o,'3 ~ + (^w{Cq3 — Cq2) + (^uCo^sft] tan/3 



^9,3 



Xfc 



(8) 



[1 + (a,Co^;| + tan/3][l + (a,Co^;| + a^Co^g + auC^Ji) tan/3] ' 

We have checked that if we multiply (7) by v^, we recover the SUSY correction to the b 
quark mass as given by Pierce et al. [25] . For degenerate squarks, we recover the expression 
for Xfc in R-ef- [5] by setting = in the denominator of (8). 

The flavour violating couplings between first and third generation down squarks are ob- 
tained by obvious substitutions. If squarks of the first and second generations are degenerate 
(as is the case in many models, including all models we consider in this paper) we can replace 
sl by without changing Xfc- Therefore, for the calculation we only have to replace 
Vts Vtd in Eq. (6), keeping everything else the same. 

Eq. (8) ignores intra-generation squark mixing. Assuming that this is significant only 
for third generation squarks, we find the modified result. 



Xfc = -— tan/3 , 



(9) 



with 



Co{mg, mi^,mij - slCo{mg, mi^.m-s^) - clCo{mg, mi^,msj + auftCo{ii, mi^.mQ 



L> = |l + agtan/3 slCo{mg,mi^,msj^) - clCo{mg,mi^,msj^) 



+ a^tan/3Co(M2,/x, ma^) 



X |l + tan l3Co{mg, m^^, m^^) + tan /3 
+auft tan /3Co (//, m^^ , m^^ ) | , 



c^Co{M2, n, mi) + slCo{M2, n, m 



t2J 



(10) 



where Sft^j = sin 6'^^^, Cft^^ = cos^^^j, and Oi^^t are the bottom and top squark mixing angles, 
respectively. 

We should note that our calculation includes only terms that are most enhanced by 
powers of tan/3. Our calculation is, therefore, valid only when tan/3 ^ 25 — 30. We will 
subsequently see that for smaller values of tan /3 the branching fraction of interest is too low 
to be of interest at the Tevatron, so that this is not a serious handicap for the purpose of our 
analysis. We remark that in the loop functions we have assumed that the chargino masses 
are well approximated by IM2I and 

Finally, the Icptonic couplings of the neutral Higgs bosons h, H and A of the MSSM 
that are needed to complete the evaluation of Higgs boson mediated leptonic decays of Bg 
or B^ are given by. 



■'HU 



[h sma + H cos a\ u H i'fsiA . 



2Mw cos /3 



2Mi 



w 



(11) 
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B. The branching ratio for Bg^a decays 

The effective Hamiltonian [15] for B^' — > i'^i^, d! = s, d, is given by 

G 

H = -^Vi^,Vtb[cwOio + cq^Qi + cq^Qs] + h.c. , 



witfi the relevant operators being given by, 



;d'LbRie, Q2 



d'Lhdld 



The coefficient cio that encapsulates the SM contribution is given by [15], 

Y{xt) 



ClO 



-4.2 



(12) 



(13) 



The function Y{xt) that appears here is defined in Ref. [15]. The SUSY contributions can be 
obtained by applying the matching condition between the amplitude given by the effective 
Hamiltonian and the one obtained using (6) and (11). We find. 



27r 

cqi = — Xfc 



mi,me /cos(/3 + a) sin a sin(/3 + a) cos a 



CQ2 



a '^^'^cos^/3sin^/3 
27r mbiTii 1 

X,FC 9 r> 2" ■ 

a cos"^ p m\ 



mt 



(14) 



Note that as — > 00, a + /3 — > | and mjj ~ m^. This then implies cq^ ~ "Cqj, an 
observation that will prove useful later. 

Finally, using the hadronic matrix elements. 



{0\h^^,d'ix)\Ba'iP))^tfB„P^e 
{Q\h^d\x)\Ba>{P)) = -ifB.- 



-iP-x 



-iP-x 



we can write the branching ratio for the decay Bd' — > d'^d^ as. 



BiBd' 



647r3 



X 



1 - 



Ami 



ml 



rrib + md' 



Amj 



m 



+ 



2mf 



mB 



-cio 



■d' 



(15) 



C. Comparision with other studies 

A complete calculation of the 1-loop chargino-induced fiavour violating MSSM Higgs 
boson couplings, along with the corresponding result for the coefficients Cq^ and Cq^ in 
(14), may be found in Ref. [8]. In comparision, we have only retained leading terms in 
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tan /3, so that our simplified calculation is valid only if tan /? 25 — 30. We have checked 
that our result for the contribution to the coefficients Cq^ ,^ from the chargino graphs indeed 
reduces to (5.13) of Ref. [8], with ^^y^ ^ ~ ^^^^ The results of Ref. [8] have been 
used for the recent phenomenological analysis [10,11] within the mSUGRA framework.^ Our 
calculations improve upon these in that we include contributions from the gluino-mediated 
flavour violating contributions (which we will see are significant) first discussed by Babu and 
Kolda [5]. We found it difficult to tell whether or not these contributions were included in 
the recent phenomenologial analysis of the mGMSB and mAMSB scenarios by Back et al. 
[14] : they do not give any formulae nor do they discuss how their calculation was performed. 

While we were preparing this paper, Ref. [12] appeared. This study includes a complete 
calculation of gluino and also analogous neutralino-induced contributions to fiavour violating 
Higgs boson couplings, but again within the framework of the MSSM with specific assump- 
tions about sfermion mass matrices. It is pointed out that it is possible to find special regions 
of MSSM parameter space where due to large cancellations between gluino and chargino con- 
tributions, neutralino diagrams (which are usually small) can no longer be neglected. Indeed, 
in this case even the ratio of the branching fractions B{Ba — > pL'^ii')/ B{Bs — > A*^A*~) may 
be different from that expected from ratios of Kobayashi-Maskawa matrix elements. 

III. B(Bs,d ^ ^+^~) IN SUSY MODELS: RESULTS 

We are now ready to proceed with the evaluation oi BiyBs^d ~^ in the various SUSY 

models discussed above. For a given point in the parameter space of these models, we use the 
program ISAJET v 7.63 [26] to evaluate the corresponding MSSM parameters that we need 
for all phenomenological analysis, including the computation of this branching ratio. This 
then allows us to assess the ranges of model parameters where signals from B{Bs (i £~^i~) 
decays may provide evidence of deviation from the SM. A clear advantage of this approach 
is that for any given model, it is also possible to compute SUSY contributions to other 
observables; e.g. B{b — > 57) or g^ — 2 that have already been constrained by experiment 
[27,28]. Moreover, constraints from the non-observation of any sparticles [29] or Higgs bosons 
[30] which yield lower limits, m^^ > 103 GeV, mg > 100 GeV, rrif, > 76 GeV, ruh > 113 GeV 
and mA > 100 GeV, on their masses can be readily incorporated.^ 

ISAJET 7.63 includes several improvements over previous versions. From our point of 
view, the most important of these is the improvement of the bottom Yukawa coupling that 
enters the computation of m^. We see from (14) and (15) that the value of plays an 
important role in the determination of the branching ratio of interest. There are regions 



^The earlier studies [7,6,9] were performed within the MSSM for particular choices of parameters. 

^Although the exact limits on sparticle masses depend somewhat on the model, and also on where 

we are in parameter space, the limits indicated here are applicable for a wide class of models. The 
lower limit of 114 GeV on the mass of the SM Higgs boson has to be translated into the limit on 
nih- Except when A is light, h is essentially the SM Higgs boson and this limit applies essentially 
without modification. 
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of parameter space where, because of cancellations, uia may be considerably smaller than 
other sparticle masses. Especially for these parameter ranges, the improved computation of 
the Yukawa coupling plays a crucial role. 

The numerical values of various meson masses, lifetimes, decay constants and Kobayashi- 
Maskawa mixing matrix elements that we use as inputs to our analysis are, 

Bg mesons: 

niB, = 5.3696 GeV, /b, = 0.250 GeV, tb, = 1.493 ps 

Bd mesons: 

ruBa = 5.2794 GeV, /b, = 0.208 GeV, tb^ = 1.548 ps 
Kobayashi-Maskawa matrix elements: 

\Vtb\ = 0.999, \Vts\ = 0.039, \Vtd\ = 0.009 



A. Minimal supergravity model (mSUGRA) 

SUSY is assumed to be broken in a hidden sector consisting of fields that interact with 
usual particles and their superparners only via gravity. SUSY breaking is communicated to 
the visible sector via gravitational interactions, and soft SUSY breaking sparticle masses and 
couplings are generated. Without further assumptions, the scalar masses can be arbitrary 
leading to flavor changing processes in conflict with experiment. 

Within the mSUGRA grand unified framework [18], it is assumed that at some high scale 
(frequently taken to be ~ Mqut) all scalar fields have a common SUSY breaking mass mo, 
all gauginos have a mass mi/2, and all soft SUSY breaking scalar trilinear couplings have a 
common value Aq. Electroweak symmetry breaking is assumed to occur radiatively. This 
fixes the magnitude superpotential parameter fi. The soft SUSY breaking bilinear Higgs 
boson mass parameter can be eliminated in favour of tan so that the model is completely 
specified by the parameter set: 

^0-1 ^1/2-1 ^0, tan/3, sign{ii). (16) 

The weak scale SUSY parameters that enter the computation of sparticle masses and cou- 
plings required for phenomenological analyses can be obtained via renormalization group 
evolution between the scale of grand unification and the weak scale. 

We use the program ISAJET to evaluate these MSSM parameters, and via these the 
various masses and mixing angles that enter the fiavour violating coupling xfc given by 
Eq. (10). The required partial width can then be computed using Eq. (15). 

In Fig. 1 we show the dependence of B{Bs — > fJ^~^fJ^~) and B{Bfi t~^t~) on a) tan/? 
for Ao = 0, and b) Aq for tan/3 = 40. We fix niQ = mi/2 = 300 GeV, and illustrate the 
results for both positive (solid) and negative (dashed) values of /x. Values of tan/3 {\Ao\) 
larger than the corresponding value denoted by squares on the curves are where a sparticle 
or Higgs boson mass (in this case, it is always m^j) falls below its experimental bound. The 
uppermost set of the dashed and solid lines corresponds to the branching fraction for the 
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decay b — >■ 57, the scale for which is given on the right hand axis. Our purpose in showing 
this figure is to understand the behaviour of the leptonic decays of and Bg, so that the 
fact that B{b — > sj) appears to be outside the allowed range should not bother the reader. 
Several features are worth noting. 

1. Since we have neglected contributions from neutralino loops [12], the branching ratio 
for Bd — > r^r~ decays is larger than that for Bs A*^/^" decays by just a constant 
factor fixed by SM parameters. 

2. Focussing for the moment on the low tan/5 end of the B{Bs fJ'^fJ' ) curves, we see 
that the solid and dashed lines lie on either side of the SM value, refiecting the fact 
that the sign of the SUSY contribution fiips with the sign of and its interference 
with the SM contribution goes from destructive for < to constructive for positive 
values of /x. Of course, this becomes irrelevant for very large values of tan/3 where the 
SUSY contribution dominates. 

3. A striking feature in both frames is the very sharp rise in the branching ratio for /i < 0. 
Naively, we expect the SUSY contribution to behave as tan^ j3 / m\7 The dashed curves 
in frame a) indeed roughly show this scaling behaviour as long as we stay away from 
the lower range of tan /? values where interference with the SM amplitude is significant. 
But this is somewhat fortituous because the dashed curves in frame b) rise much faster 
than would be naively expected. This is because it is also very important to take 
into account differences in other MSSM parameters, most notably At and fi that enter 
via Qu and in (3). In general, the dashed curves are steeper than the solid curves 
because, for n < 0, rriA decreases sharply as tan/3 increases: there is no corresponding 
decrease of positive values of /i. 

4. It is interesting to sec that both B{Bs /^"""/x") and B{Bd t^t^) can be potentially 
very large (and even exceed the current upper limit) for values of parameters where 
there are no direct signals for SUSY or Higgs bosons. This region of parameters is also 
excluded by the experimental value of the inclusive branching ratio B(b ^ 57) [27]. 
It should be understood though that our purpose in showing this is pedagogical. 

5. We have also examined the case with rriQ — 1 TeV and other parameters as before, 
except that the allowed range of Aq is much larger. For positive values of ^ and Aq = 0, 
B[Bs Ai"*"/^") varies between (4 — 100) x 10~^. For < 0, this branching fraction 
is close to or smaller than the SM prediction, except when 46 < tan /3 ^ 48 when the 
branching fraction shoots up by over two orders of magnitide. For tan /3 = 40 as in 
Fig. lb and A* > 0, the corresponding dependence on is qualitatively similar to that 
the figure. For negative values of /j, the branching fraction is close to its SM value for 
\Ao\ ^ 1.2mo, but rapidly shoots up to beyond 10^^ when \Ao\ becomes large. As in 
the mo = 300 GeV case, the rapid rise is driven by the rapid drop in m^. 



^Recall that the contribution from h exchange is very small as long as /i is a SM-like Higgs, and 
further, that itih ~ tua in the same limit. 
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Having used this pedagogical illustration to obtain an idea of how and why the branching 
fractions vary with parameters, we now turn to an exploration of these branching ratios 
in mSUGRA parameter space. Because the sparticle spectrum is most sensitive to mo 
and mi/2, the mo — mi/2 plane provides a convenient way of displaying the results, as 
illustrated in Fig. 2 for ^o = and a) /j, > 0, tan/3 = 50 and b) /j, < 0, tan/3 = 45. The 
dark-shaded regions are excluded by theoretical constraints: charge-breaking minima or 
lack of electroweak symmetry breaking. In the slant-hatched region, Zi is not the lightest 
supersymmetric particle (LSP). The region covered by open circles is excluded by lower 
limits on sparticle masses or on m^. Below the dot-dashed hne labelled h, ruh is smaller 
than 113 GeV: we show this separately because this limit is modified if m^ < 150 GeV. Our 
main results are the contours of B{Bs — > IJ''^IJ'~) (solid) and B{Bd — > r+r~) (dashed) within 
the mSUGRA framework. The contours are labelled by the values of the corresponding 
branching fractions. In frame a), the innermost dashed contour corresponds to a branching 
fraction of 10~^. Prom frame a) we see that even in the region allowed by all these constraints 
B(Bs — > lJi^lJi~) {B{Bd — > r'^r~)) may be as large as 10~^ (10~^), while from frame b) the 
corresponding branching fraction may be considerably larger. Also shown are contours of 
fixed branching fraction for the decay 6 — > S7 computed using the calculation of Ref. [31] 
with improvements described in Ref. [17]. Below the contours labelled 6 — > 57 in frame a) 
(frame b)) the branching fraction is smaller (larger) than 2xl0~^ (5 x 10"'^), and appear to 
be comfortably outside our assessment 2.16 x 10~^ < S(6 — > S7) < 4.34 x 10~^, [17] for the 
allowed experimental range [27] of this branching fraction. We caution the reader that the 
SUSY contribution, for large values of tan f3, may have considerable theoretical uncertainty, 
so that this "constraint" should be interpreted with some care.^ 

The dotted curves arc contours of B{Bs — * /^^/^^) = 10~^, obtained by calculating the 
branching ratio by retaining just the chargino-loop in evaluating the SUSY contribution 
to the decay i.e. if (ig in (3) is set to zero. We remind the reader that it is just this 
contribution has been included in recent analyses [10,11]. We have checked that the dotted 
contour in frame a) is in good agreement with the corresponding contour in Fig. 3 of Ref. [10]. 
This provides quantitative confirmation of the validity of our approximations, relative to the 
complete calculation of Ref. [8] . More importantly, the difference between the dotted contour 
and the corresponding solid contour highlights the importance of retaining the effect of gluino 
loops whose contribution appears to interfere destructively with the chargino contribution. 
Thus, earlier conclusions [10,11] about the SUSY reach of Tevatron experiments may be 
over-optimistic. 



^It is also worth emphasizing that unlike constraints from direct searches, constraints from B{b 
S7) are very sensitive to details of the model. For instance, small amount of flavour mixing (from 
some unknown physics) in the squark sector could lead to large differences in the predictions for 
b — > 57, and for that matter, Bg' — > £^i~ decays. For this reason, we urge our experimental 
colleagues to view theorists' assessment of "excluded regions" (especially when these are excluded 
due to SUSY loop effects as opposed to direct searches) including those in this paper in the 
proper perspective. It is logically possible that small deviations from the defining assumptions of a 
particular framework such as mSUGRA may permit much larger signals without being in conflict 
with current constraints. 
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To see how the sensitivity of these experiments varies with tan /3, we show in Fig. 3 
contours of B{Bs A^^/^ ) = 10^'' (sohd hnes) and B{B,i r+r") = 10^^ (dashed 
hnes) for the values of tan (3 that label the contours. We take Aq = and show results 
for a) fj, > and b) /i < 0. The hatched region corresponds to tan/3 = 35. The lines in 
frame a) terminate at the corresponding boundaries of the theoretically forbidden regions. 
As expected, the region over which the branching fraction may be probed at the Tevatron 
is very sensitive to tan/3. Even for tan/3 = 35, this channel appears to yield a better SUSY 
sensitivity than direct searches [32] for an integrated luminosity of ~ 2 fb~^. 

Prom Fig. 2 and Fig. 3 we note the following: 

• The branching fraction under study is significantly larger for negative values of fi. We 
have traced this to the structure of the denominator of xfc in Eq. (10). Changing the 
sign of II changes the sign of all the a^, so that a suppression for positive /i changes to 
an enhancement for // < 0. Notice that while our denominator factor reduces to that 
in Ref. [5] in the appropriate limit, it differs from that in Ref. [10] who do not have 
the ag term in the second factor.^ Quantitatively, this leads to differences of 0(10%) 
and do not affect the qualitative conclusions. 

• Although /X < is generally thought to be disfavoured by the determination of the 
muon anomolous magnetic moment by the E821 experiment [28], we advise caution 
in this regard: in contrast to conventional wisdom [33] a conservative estimate [34] 
of the theoretical error suggests that there is a region allowed [17] by this constraint, 
though perhaps in conflict with B{b — > 57), where Bs — > decays may provide 
the flrst hint of new physics if tan /3 42 or so. This region would expand as Tevatron 
experiments accummulated more data. 

• Prom these flgures, we see that a sensitivity of 10~^ for B(Bd — > r+r") would roughly 

yield the same reach as an order of magnitude better sensitivity that is expected to be 
attained at the Tevatron with an integrated luminosity of ~ 2 fb^^. We do not know 
whether experiments at B factories will be able to attain this sensitivity. 

Up to now, in our scan of mSUGRA parameter space, we have only considered models 
with Aq = 0. To understand just how large the branching ratio for Bg fi^fi^ can become 
for other values of Aq, we have scanned mSUGRA models with mo between 100 GeV and 
-\/50mi/2 and — 3mo < Aq < Suiq. In Fig. 4, we show B{Bs — > A* "*"//") as a function of mi/2, 
for a) tan/3 = 55, > 0, b) tan/3 = 45, /i < 0, and c) tan/3 = 40, // < 0. Models are 
accepted only if they satisfy the theoretical and the experimental constraints from direct 
searches for sparticles and Higgs bosons. For each model, we put a cross (dot) if i3(6 — >■ 57) 
lies within (outside) (2 — 5) x 10""^. We note the following: 

1. Although we have shown just B{Bs — » ^'^ in this figure, the range of B[B^ — > 
T+T"), which differs only by a constant factor can easily be estimated. 



^In Ref. [8], which is a strict diagramatic calculation, does not have this denominator. Dedes et 
al. who use the results of Ref. [8] , include a denominator correction in their formulae. 
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2. We see that for > 0, in frame a) the branching fraction is typically smaller than 
10~^, and never exceeds ~ 3 x 10^^ for models loosely satisfying the b ^ constraint. 
But for this constraint, it can be as high as 10~^. 

3. Turning to frame b) we see that for negative values of fi the branching fraction can be as 
large as 10~^, but most of these points are excluded by the experimental measurement 
of 6 — > S7 because the predicted branching fraction is too large. For large values of 
mi/2, however, B{Bs —>■ may even exceed 10~^ for parameter values allowed by 
b ^ constraints. For the slightly smaller value of tan/3 shown in frame c), we see 
that these large values of B{Bs A*"^A*~) are possible over a wide range of mi/2. 

4. In all frames we see that although the branching ratio may be much larger than the 
SM value for some set of SUSY parameters, it never falls below about half the SM 
value. Experiments at the LHCb should be sensitive (at the 3o" level) down to this 
branching fraction [35], so that a non-observation of this decay in these experiments 
would signal new physics other than the scenarios considered here.^° 

5. A striking feature of the figure is the "gap" for small mi/2 values in frame a). It appears 
that for positive ^ and very large values of tan/5 the branching ratio in mSUGRA 
is unlikely to be at the SM value. Of course, when mi/2 (and hence m^) is very 
large, the branching ratio tends to the SM value. The gap is absent for // < 0. The 
existence of this gap can be qualitatively understood by focussing on the chargino 
contribution (proportional to the term) and recognizing that models with positive 
/i allow both signs of Aj, while models with negative almost always have < 0. As 
a result, the SUSY contribution may interfere constructively or destructively with the 
SM contribution when // > 0, while for negative values of // the interference is almost 
always destructive. Of course, the interference is important only when the SUSY and 
SM amplitudes have comparable magnitudes. For the large values of tan/3 in this 
figure, this happens only if is large (otherwise the SUSY amplitude is much larger 
than the SM one). In this case, the coefficients cq^ and cq^ in Eq. (15) have the same 
magnitide.^^ In units where the SM contribution in the square parenthesis of Eq. (15) 
is 1, if we write the cq^ contribution as x = —CQ^mBj^ciom/j^, the partial width is 
determined by F = + (1 + a;)^. This has a minimum of 1/2, explaining why the 
branching fraction does not fall below about half its SM value. Moreover, as long as x 
is positive and not small (as is the case for small values of mi/2) F significantly exceeds 
the SM value, li x < 0, F becomes large only if \x\ is very large, thereby accounting 
for the gap for positive values of /i. For negative values of /i, there is no "positive x 
branch" , and hence, no gap. 



Admittedly we have not scanned all parameter space. However, if tan/3 is small the SUSY 
contribution is reduced, and we would expect the branching fraction to be closer to its SM value. 

^^This really follows because of the properties of the MSSM Higgs sector, and so should be true 
in other models also. 
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6. Although wc have not shown this, the scatter plot for tan (3 — 50, n > is very similar 
to frame a), except that the "gap" is somewhat less pronounced. In other words, the 
difference between the distribution of dots and crosses in frames b) and c) is absent. 

B. Minimal gauge-mediated SUSY breaking model (mGMSB) 

In these models, SUSY breaking is again assumed to occur in a hidden sector which 
somehow couples to a set of messenger particles that couple not only to this hidden sector, 
but also have SM gauge interactions. The coupling to the hidden sector induces SUSY 
breaking in the messenger sector, which is then conveyed to superpartners of usual particles 
via SM gauge interactions. Messenger particles are assumed to occur in complete vector 
representations of SU{5) with quantum numbers of SU{2) doublets of quarks and leptons. 
We assume ns < 4 since otherwise gauge couplings do not remain perturbative up to the 
scale of grand unification. The messenger sector mass scale is characterized by M. The 
soft SUSY breaking masses for the SUSY partners of SM particles are thus proportional 
to the strength of their corresponding gauge interactions, so that squarks are heavier than 
sleptons. Gaugino masses satisfy the usual "grand unification" mass relations, though for 
very different reasons. Within the minimal version of this framework, the couplings and 
masses of the sparticles in the observable sector are determined (at the messenger scale M) 
by the parameter set, 

A, M,n5, tan /3,sign{n),C grav (17) 

A sets the scale of sparticle masses and is the most important of these parameters. The 
parameter Cgrav > 1 and enters only into the partial width for sparticle decays to the 
gravitino and is irrelevant for our analysis. The model predictions for soft-SUSY breaking 
parameters at the scale M are evolved to the weak scale using ISAJET and used to compute 
B{Bg, i+i-) as before. 

The novel feature of the GMSB framework is that SUSY breaking can be a relatively low 
energy phenomenon if the messenger scale is small. It is in this case that the gravitino can 
be an ultra-fight LSP. In such a scenario heavy sparticles cascade decay to the next lightest 
supersymmetric particle (NLSP) which then decays to the gravitino and an ordinary particle. 
The NLSP may be the lightest neutralino or the stau.^^ SUSY signatures at the Tevatron 
[36] are sensitive to the nature of the NLSP. In our considerations we will focus on models 
with relatively low messenger scales for which the collider phenomenology differs the most 
fi-om mSUGRA. 

An important difference between this framework and the mSUGRA model is that A- 
parameters (at the scale M) which are generated only at two loops are much smaller than 
scalar and gaugino masses. The weak scale parameter Af that is obtained by rcnormaliza- 
tion group evolution from the messenger scale is typically smaller than in mSUGRA if the 
messenger scale is not large. We thus expect that that sets the scale of the chargino 



For low tan/3 values not of interest to us, the sleptons of different generations are essentially 
degenerate leading to the so-called co-NLSP scenario [36]. 
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contributions in (3) is significantly smaller than in the mSUGRA framework. We also ex- 
pect that the ghiino contribution to the loop decay of Bg/ would also be smaller than in 
mSUGRA because there is not enough room to run and induce large mass splitting between 
the different Qls by renormalization group evolution if the messenger scale is low. There- 
fore, we typically expect lower values of B{Bqi — > i'^£~) relative to corresponding results for 
mSUGRA in these scenarios. This branching fraction would, for the same reasons, increase 
logarithmically with increasing messenger scale. 

Our results for the branching fraction within the GMSB framework are shown in the A — 
tan l3 plane in Fig. 5. Here, we fix M = 3A, — 1 and illustrate the results for both signs of 
ji. As before, the dark-shaded region is excluded by the theoretical considerations on EWSB 
discussed previously. The bulge on the left comes from m|j < 0. while in the horizontal 
strip at high tan j3 that extends to large A values, m\ < 0. In the slant-shaded region 
rufj^ < 76 GeV, while in the slivers covered by triangles along the boundary of the horizontal 
dark-shaded region, tha < 100 GeV. Finally, values of A up to A ~ 53 TeV (covered by 
open circles) are excluded because mg^ < 100 GeV, and possibly also by constraints on the 
chargino mass or decay properties [37]. To the left of the dot-dashed contour labelled 
h, rrih < 113 GeV. The branching fraction for the decay 6 57 is between (2 — 5) x 10~^ 
over the entire plane in frame a), while it is in this range to the right of the corresponding 
contour in frame b). Contours of B{Bs fJ^~^l^~) and B^B^ — > r+r") are shown as solid and 
dashed hues, respectively. These are labelled by the value of the corresponding branching 
fraction. We see that over almost the entire plane B{Bs —>■ {B{Bd t~^t~)) is 

smaller than 10~^ (10~^) confirming our earlier expectation that the branching fractions 
will be smaller than in mSUGRA. It would be difficult to probe these decays at this level 
at existing facihtics. It is only for the largest values of tan P where we approach the region 
where dives to very small values that these branching fractions might be accessible: 
however, of this region in frame b) is "excluded'' by the experimental value of B{b — > 57). 
In contrast, direct SUSY searches should be sensitive to A values of 118-145 TcV depending 
on the integrated luminosity that is accumulated [36,38]. Direct searches at the LHC will 
easily probe the entire plane [39]. 

Up to now, we have assumed — 1. Since scalar (gaugino) masses scale with (77,5) 
we would, expect large sensitivity to if we show the results in terms of A. To avoid a 
proliferation of figures, in Fig. 6 we show the branching fractions versus rrig for 71,5 = 1 — 4. 
We have fixed M = 3A and tan/? = 40 and illustrated the results for both signs of fi. 
The curves are terminated at values of rrig that violate any of the constraints discussed 
previously. For very large values of rrig the branching fractions approach the SM values 
because tends to be large. For smaller values of rrig the branching fractions (for a given 
value of rrig) do depend on n^. Nevertheless, it seems that the range over which these vary is 
insensitive to the choice of n^, so that our general conclusions drawn from parameter scan 
for ns = 1 should, broadly speaking, remain unaltered. 

C. Minimal Anomaly-Mediated SUSY Breaking Model (AMSB) 

Within the supergravity framework, sparticle masses always receive loop contributions 
originating in the super- Weyl anomaly when SUSY is broken. These loop contributions are 
generally much smaller than tree level masses. There are, however, classes of models where 
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these loop contributions may dominate [20]. These include models where there are no SM 
gauge singlet superfields that can acquire a Planck scale vev and the usual supergravity 
contribution to gaugino masses is suppressed by an additional factor ^^']^/^^ relative to 
1713= MgugY / ^P-i higher dimensional models where the coupling between the observable 
and hidden sectors is strongly suppressed. 

The anomaly-mediated SUSY breaking contributions to gaugino masses are proportional 
to the corresponding gauge group /9-function, and so are non-universal. Likewise, scalar 
masses and trilincar terms arc given in terms of gauge group and Yukawa interaction beta 
functions. As a result sparticles with the same gauge quantum numbers have a common 
mass so that flavour changing effects from this sector are naturally small. Slepton squared 
masses, however, turn out to be negative (tachyonic). A "minimal" fix, that does not upset 
the resolution of the flavour problem is to assume an additional contribution rn^ for all 
scalars. Since the magnitude of ^ is flxed by the observed value of M^, the parameter space 
of the model then consists of 

"^0) "^3/2, tan/3 and sign{^). (18) 

Weak scale SUSY parameters can now be computed via renormalization group evolution, 
and B{Bqi — > can again be calculated using (15). 

For small values of mo, the scale of sparticle masses is set by m3/2 which cannot be 
much smaller than about 35 TeV. In Fig. 7 we show the branching ratio for Bs fi^fi~ 
and Bd — >■ t~^t~ decays versus tan /3, with 777.3/2 = 40 TeV and a) mo = 400 GeV, and 
&) 777o = 1 TeV. The sohd (dashed) curves correspond to positive (negative) values of fi. As 
before, the squares mark the experimentally allowed upper limit on tan (3. Also, shown by 
the right hand scale in the flgure, is the branching fraction for the decay b — > 57. 

We see from Fig. 7 that for positive values of fi the B{Bqi £^i~) is close to its SM value, 
and relatively insensitive to the value of tan /3. For negative values of /7, these branching 
fractions can indeed become large, again when niA dives to low values. However, for the 
range of tan/3 where B{b — > 57) < 5 x 10~^, the branching fraction for Bg — > A^''"/^" decay 
will be difficult to detect in Tevatron experiments. We have also examined the dependence 
of these branching fractions on 7/73/2. For positive values of fi, and tan /3 = 40 — 60, they are 
once again always close to the corresponding SM value, independent of 7773/2. For negative 
values of /i, B{Bs /7+/7^) increases rapidly with 7773/2, and even for tan/3 = 40 and 
777o = 1 TeV can be as large as 10~^; however, for the range of 7773/2 where S(6 — > 57) is not 
too large, the branching fraction can be as large as several times 10~^ and may be accessible 
at the Tevatron, but would require considerable integrated luminosity. 

In Fig. 8, we show contours for B{Bq/ — > £+£~) in the 7770 — 7773/2 plane, where we have 
fixed tan /3 = 42. We show results only for < since for positive values of /j, the branching 
fraction appears to be close to the SM one even for very large values of tan /3. We remind 
the reader that in the AMSB framework negative fi is also favoured by the result of the 
E821 experiment [28]. The dark-shaded region is excluded by theoretical constraints and 
in the open circle region itl^^ < 100 GeV. Along the line of squares and triangles running 
diagonally across the figure, nih or rriA, respectively fall below their experimental bounds. 
Finally, except in the narrow region covered with dots that follows the boundary of the 
upper shaded region where B[b ^ S7) is too large, the branching fraction for the decay 
6 — > 57 is in its "allowed range" of (2 — 5) x 10~^. The three sohd contours correspond 
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to B{Bs — >■ IJ-'^IJ-~) = 10~^, 5 X 10~^ and 10~^, while the dashed contours correspond to 
B{B(i r+r") = 10^'^, 10^^ and 10^"''. We see that the branching fractions within the 
AMSB framework is small over most of the parameter plane. The decay Bg — > fJT 
may be probed in Tevatron experiments only over the limited range where niA is diving 
to relatively low values. With a large integrated luminosity, there is a limited parameter 
range (consistent with other collider constraints) where this decay may be the harbinger of 
new physics at the Tevatron. With just 10 fb^^ of integrated luminosity direct searches for 
sparticles at the LHC should be sensitive [40] to 777.3/2 ~ 70 — 90 TeV, depending on the 
value of 7770- 

For larger values of tan (5 and < 0, the situation is qualitatively similar except that 
the region excluded by the theoretical constraints expands leaving a yet narrower wedge of 
"allowed parameters" . 



D. SU{5) models with non-universal gaugino masses 

Since supergravity is not a renormalizable theory, there is no reason to suppose that the 
gauge kinetic function fab — Sab- Indeed, if the gauge kinetic function develops a SUSY 
breaking vev that also breaks the GUT symmetry, non- universal gaugino masses result. 
Expanding the gauge kinetic function as fat = Sab + ^ab/Mpianck + • • •, where the fields ^ab 
transform as left handed chiral superfields under supersymmetry transformations, and as 
the symmetric product of two adjoints under gauge transformations, we parametrize the 
lowest order contribution to gaugino masses by, 

jCd [ d^eW'^W'-p^ + h.c. D i^il^A^A^ + . . . , (19) 

J -Mpianck -Mpianck 

where W"' is the supcrfield that contains the gaugino field A", and F$ is the auxiliary field 
component of $ that acquires a SUSY breaking vev. In principle, the chiral supcrfield $ 
which communicates supersymmetry breaking to the gaugino fields can lie in any represen- 
tation contained in the symmetric product of two adjoints, and so can lead to gaugino mass 
terms that break the underlying gauge symmetry. We require, of course, that SM gauge 
symmetry is preserved. 

In the context of SU{5) grand unification, F$ would most generally be a superposition 
of irreducible representations which appears in the symmetric product of two 24s, 

(24x24)sy^^etric = 1 e 24 ® 75 e 200 , (20) 

where only 1 yields universal gaugino masses as in the mSUGRA model. The relations 
amongst the various GUT scale gaugino masses have been worked out e.g. in Ref. [21], and 
are listed in Tabic I along with the approximate masses after RGE evolution to Q ~ Mz- 
Motivated by the measured values of the gauge couplings at LEP, we assume that the 
vev of the SUSY-preserving scalar component of $ is neglible. Each of these three non- 
singlet models is as predictive as the canonical singlet case, and all are compatible with the 
unification of gauge couplings. Although superpositions are possible, for definiteness we only 
consider the predicitive subset of scenarios where F$ transforms as one of the irreducible 
representations of SU{5). The model parameters may then be chosen to be. 
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mo, M3, Aq, tan P and sign{fj,) , 



(21) 



where M° is the SU (i) gaugino mass at scale Q = Mqut'- -^2 ^^'^ ^1 then be calculated 
in terms of using to Table I. The sparticle masses and mixing angles can then be 
computed using the non-universal SUGRA option in ISA JET. An illustrative example of 
the spectrum may be found in Ref. [41] where this framework is reviewed. We see from 
Table 1 that the pattern of gaugino masses at the weak scale differs quite significantly from 
mSUGRA expectation. This suggests that the relative contribution of the chargino and 
gluino-mediated contributions to Bg tJ'^ tJ'~ decays, and hence the branching fraction may 
differ significantly from their values in the mSUGRA framework. 

To illustrate this, we show B{Bs f^^fJ' ) as a function of the input GUT scale gluino 
mass parameter in Fig. 9 for the various models in Table 1. We fix mo = 400 GeV, Aq = 0, 
tan /5 = 40 and show our results for both signs of /i. The dotted line labelled 1 DR shows 
the branching fraction within the mSUGRA framework where $ transforms as a singlet of 
SU (5), while the solid and dashed lines show the result for the cases where l» transforms as 
a 24 or 200 dimensional representation, respectively. The 75 case is theoretically excluded 
for our choice of parameters. Indeed, we sec that the branching fraction is sensitive to the 
underlying pattern of non- universality. Especially striking is the dip at M3 ~ 200 GeV 
in frame b). We have checked that this is due to an almost complete cancellation between 
the chargino and gluino loops, which as far as we can ascertain is completely accidental. 
Individually, either the chargino or gluino loop contribution would be much larger than the 
SM one (m^i is only 175 GeV), but the total SUSY amplitude is comparable to the SM one 
and the two interfere destructively to yield a branching fraction which is close to its minimum 
value (recall our discussion at the end of Sec. Ill A). This case highhghts the importance of 
including the gluino loop, without which the branching fraction would be almost an order 
of magnitude bigger. We should also mention that contributions from neutralino loops that 
have been neglected in our analysis may now be significant, in which case B[Bd t^t~) 
would not have the canonical ratio with B[Bs — >■ [12]. 

To see how large B{Bs — > Z^^/^^) could be in these scenarios, we scanned over the 
parameter space of these models. We allow mo to vary between 100 GeV and \/50M^ and 
—3mo < Aq < 3Ao. In Fig. 10, where we take tan/3 = 55, we show results for > 
for the three cases of non-universality discussed above. The corresponding result for the 
mSUGRA case was shown in Fig. 4a. Analogous results for < and tan (3 = 45 are shown 
in Fig. 11 with the corresponding result for the mSUGRA case now in Fig. 46. Again, points 
that satisfy other experimental constraints-*^^ but where S(6 — > 57) falls within (outside) the 
range (2 — 5) x 10~^ are denoted by a cross (dot). The following features of these figures are 
worth noting: 



^^Values of smaller than 200 GeV (300 GeV) are "excluded" for the 200 (24) model since 
these yield mh < 113 GeV. 

^^In the case of the 75 and 200 models the mass gap between the chargino and the LSP is small, 
and LEP constraints on cited here may be too stringent; see, however, the discussion in Ref. 
[42]. 



19 



1. For the 24 model, the branching fraction is always smaller than about 2 x 10~^ at 
least for the tan/? values shown, and would be difficult to probe at the Tevatron. 
However, in this scenario that Tevatron experiments should see signals in the 
jets + ie+ ^r, and possibly, jets + Z° ^ £+£-+ channels if < f 75 GeV and 
mo ^ 400 — 500 GeV [42] . For n < [43] favored by the measured value of the muon 
anomalous magnetic moment, the branching fraction tends to be within a factor ~ 2 
of its SM value. 

2. The branching fraction for Bg — > /^"^A*" decays can be large for the 75 model, regardless 
of the sign of /i. For positive values of n favoured by the E821 experiment [28], 
B{Bs — > fJ'^fJ'^) can be as large as 10~^, which is just a little over a factor of 2 from its 
current upper limit. The gap in Fig. 106 is presumably for the same reason as in the 
mSUGRA case. It is, however, more pronounced because small values of \Ao\ are not 
allowed. For the 75 model, direct SUSY searches at the Tevatron can give observable 
signals only in the channel, and that too only if Mg and mo are rather small. In 
this case, "my^^ , m^^ and m^^ are all rather close in mass, so that decays of Wi and Z2 
typically lead to very small visible energy [42] . 

3. For the 200 model, B{Bs /^"'"A*") can exceed 10~^, even for very large values of 
Mg for both signs of fi and parameter values consistent with the measured value of 
B{b — > sj). As in the 75 case, direct SUSY signals at the Tevatron, for the most part, 
will be restricted to the channel [42] . 

IV. SUMMARY 

If tan f3 is large, substantial flavour violating couplings of neutral MSSM Higgs bosons to 
down type quarks are induced at the 1-loop level via diagrams with squarks and charginos or 
squarks and gluinos in the loop. These induced interactions, in turn, lead to new Higgs boson 
mediated SUSY contributions to the amphtude for the decays Bg yW^A*" and Bd — > r+r". 
Depending on model parameters, the branching fraction for the decays may be several orders 
of magnitude larger than its corresponding SM expectation. Moreover, while the SUSY 
contribution to the decay decouples when m^ becomes large (so that the Higgs sector of the 
MSSM reduces to that of the SM with a light Higgs boson) , it does not decouple for heavy 
sparticles. 

The CDF experiment has already estabhshed that B{Bs — > //"^A* ) ^ 2.6 x 10~^, and with 
the data sample of ~ 2 fb~^ that is expected to be collected after 2-3 years of Main Injector 
operation, will be sensitive to branching fractions smaller than about 10^''. The sensitivity 
will be even greater as Tevatron experiments approach their goal of ~ 15 fb~^. This then 
opens up the possibility of discovering new physics at the Tevatron even if sparticles are 
heavy and their signals from direct production below the level of observability. 

These considerations led us to examine predictions for B{Bs fi~^^~). This is, however, 
not possible within the context of the generic MSSM since SUSY amplitudes are sensitive to 
the a priori unknown flavour structure of squark mass matrices and trilinear couplings. In 
other words general predictions (as a function of sparticle masses) are not possible, and one 
has to resort to specific models. In this paper, we have examined how this branching ratio 
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within the framework of several popular models of SUSY breaking: the mSUGRA model, the 
minimal gauge mediated SUSY breaking model, and the anomaly-mediated SUSY breaking 
model. We also examined the range of this branching ratio in supergravity SU{5) models 
with non- universal gaugino masses. Specifically, we studied how the B{Bs Ai^yU^) varies 
with parameters, and delineated regions of parameter space where it might be observable 
at the Tevatron. Without making any representation about the sensitivity of BELLE and 
BABAR, we have also examined the branching fraction for the related decay Bd t^t~. 
Since the coupling of the Higgs to SM fermions is proportional to the fermion mass, the 
partial width for this decay is enhanced (relative to that for Bg — > H^H^) by a factor 
{mr/rrinY, but reduced by a factor (|T4(i|/|^ts|)^ that originates in the flavour violating 
Higgs boson vertex. Our rule of thumb is that these experiments will be competitive with 
Tevatron experiments if their sensitivity to B{Bd r^r^) is no more than an order of 
magnitude worse than the Tevatron sensitivity to B{Bs — > ii'^ 

Our main results are summarized in Fig. 2-4 for mSUGRA, in Figs. 5 and Fig. 6 for 
the GMSB model, in Fig. 8 for the AMSB model, and in Fig. 10 and Fig. 11 for SU{f,) 
models with non-universal GUT scale gaugino masses. We will not elaborate on the details 
here. Generally speaking, in the GMSB and AMSB frameworks, the branching fraction for 
this decay is significantly smaller than in the mSUGRA model. Within the GMSB and 
AMSB scenarios, the reach of Tevatron experiments via this decay will be relatively limited. 
A branching fraction close to 10~^ is possible only for regions of parameters where rriA 
tends to be small due to accidental cancellations. In supergravity models with non-universal 
gaugino masses, the size of the signal is sensitive to the gaugino mass pattern. Indeed, in 
the 75 and 200 models the signal can be very large for values of parameters consistent with 
other constraints. In all these models there are varying regions of parameters that are not 
in conflict with observations, where Bg — > A*"^//" should be observable at the Tevatron, and 
where there are no direct sparticle or Higgs boson signals. 

The decays Bs fi~^fi~ and Bd — * r+r" may also be incisive probes of SO (10) SUSY 
models with Yukawa coupling unification since these require tan j3 to be large. These models 
have recently received considerable attention [46], especially in light of the interpretation 
of the atmospheric neutrino data [47] as neutrino oscillations originating in a non-trivial 
flavour structure of the neutrino mass matrix. Since Yukawa coupling unification is possible 
only in very special regions of parameter space that require extensive scanning to find, we 
defer the analysis of B{Bqi £^£^) within this framework to a dedicated study of the 
phenomenology of these models that is currently in progress. 



^^We did not separately examine the gaugino- mediated SUSY breaking framework [44]. This 
is a model based on extra dimensions where usual matter and SUSY breaking fields reside on 
different spatially separated branes, while gauge fields live in the bulk. As a result, gauginos 
which directly "feel" the SUSY breaking develop masses at the compactification scale Mc, while 
matter scalars do not. Renormalization effects then induce scalar masses and A-parameters at 
Q = Mqut < Mc- Although some GUT scale non-universality is induced because of Yukawa 
couplings, we expect that results in this framework would be qualitatively similar to those in 
mSUGRA where mo,— ^ 0.5mi/2- For a discussion of the collider phenomenology of these 
models, see Ref. [45]. 
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Before concluding, we again emphasize that the branching fractions that we have exam- 
ined are sensitive to details of the model in that small deviations in some of the assumptions 
(about sfermion flavour structure at the high scale) that define the framework could result 
in considerable changes in the answer. Our point here is that we should use various "model 
predictions" only for guidance, or as benchmarks, but keep open the possibility that the real 
world may be somewhat different. By the same token, accurate determination of flavour 
violating processes serves as a sensitive probe of flavour structure at very high scales. The 
observation of processes such as Bg A*"*"/^" at the Tevatron or — > r+r" at i?-factories 
is important not only because it would herald new physics, but also because they may serve 
as probes of flavour physics at scales not directly accessible to experiment. 



ACKNOWLEDGMENTS 

It is a pleasure to thank C. Kolda for a conversation about Ref. [5] that got us started 
on this study. We thank H. Dreiner, U. Nierste and especially A. Dedes for correspondence 
and detailed comparisons of their results against ours. We are grateful to H. Baer and 
J. Ferrandis for their help in improvements in the calculation of vriA using ISAJET. But 
for their intervention, this paper may have been released considerably earlier. We thank 
E. Polycarpo for valuable information about the sensitivity of LHC experiments to Bg — >■ 
yU"*"/!" decays, and for providing us with Ref. [35]. This research was supported in part by the 
U.S. Department of Energy under contracts number DE-FG02-97ER41022 and by Fundagao 
de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP). 



22 



REFERENCES 



[1] E. Witten, Nucl. Phys. B188, 513 (1981); S. Dimopoulos and H. Georgi, Nucl. Phys. 
B193, 150 (1981); N. Sakai, Z. Phys. Cll, 153 (1981); R. Kaul, Phys. Lett. B109, 19 
(1982). 

[2] For recent reviews, see e.g. S. Martin, in Perspectives on Super symmetry^ edited by 

G. Kane (World Scientific), hep-ph/9709356 (1997); M. Drees, hep-ph/96 11409 (1996); 
J. Bagger, hcp-ph/9604232 (1996); X. Tata, Proc. IX J. Swieca Summer School, 
J. Barata, A. Malbousson and S. Novaes, Eds., hep-ph/9706307 (1997); S. Dawson, 
Proc. TASI 97, J. Bagger, Ed., hep-ph/9712464 (1997). 

[3] R. Hempfiing, Phys. Rev. D49, 6168 (1994); L. Hall, R. Rattazzi and U. Sarid, Phys. 

Rev. D50, 7048 (1994); M. Carena et al. Nucl. Phys. B426, 269 (1994); T. Blazek, 

S. Raby and S. Pokorski, Phys. Rev. D52, 4151 (1995); R. Rattazzi and U. Sarid, Phys. 

Rev. D53, 1553 (1996). 
[4] S. Rai Choudhury and N. Gaur, Phys. Lett. B451, 86 (1998). 
[5] K. Babu and C. Kolda, Phys. Rev. Lett. 84, 228 (2000). 
[6] P.S. Chankowski and L. Slawianowka, Phys. Rev. D 63, 054012 (2001). 
[7] C. S. Huang, W. Liao, Q.-S. Yan and S.-H. Zhu, Phys. Rev. D 63, 114021 (2001) and 

D 64, 059902 (2001) (E). 
[8] C. Bobcth, T. Ewerth, F. Kriiger and J. Urban, Phys. Rev. D 64, 074014 (2001). 
[9] G. Isidori and A. Retico, J. High Energy Phys. 11, 001 (2001). 
[10] A. Dedes, H. Dreiner and U. Nierste, Phys. Rev. Lett. 87, 251804 (2001); A. Dedes, 

H. Dreiner, U. Nierste and P. Richardson, hcp-ph/0207026 (2002). 

[11] R. Arnowitt, B. Dutta, T. Kamon and M. Tanaka, Phys. Lett. B538,121 (2002). 

[12] C. Bobeth, T. Ewerth, F. Kriiger and J. Urban, hep-ph/0204225 (2002). 

[13] A. J. Buras et al, hcp-ph/0207241 (2002). 

[14] S. Back, P. Ko and W. Y. Song, hep-ph/0205259 (2002). 

[15] G. Buchalla, A. J. Buras and M. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 
[16] F. Abe et al, (CDF Collaboration), Phys. Rev. D 57, 3811 (1998). 

[17] H. Baer et al, J. High Energy Phys. 07, 050 (2002). 

[18] A. Chamseddine, R. Arnowitt and P. Nath, Phys. Rev. Lett. 49, 970 (1982); R. Barbieri, 
S. Ferrara and C. Savoy, Phys. Lett. B119, 343 (1982); L. J. Hall, J. Lykken and 
S. Weinberg, Phys. Rev. D 27, 2359 (1983); for a review, see H. P. Nilles, Phys. Rep. 
110, 1 (1984). 

[19] M. Dine and A. E. Nelson, Phys. Rev. D 48, 1277 (1993); M. Dine, A. E. Nelson and 
Y. Shirman, Phys. Rev. D 51, 1362 (1995); M. Dine, A. E. Nelson, Y. Nir and Y. Shir- 
man, Phys. Rev. D 53, 2658 (1996); for a review, see G. F. Giudice and R. Rattazzi, 
Phys. Rep. 322, 419 (1999). 

[20] L. Randall and R. Sundrum, Nucl. Phys. B557, 79 (1999); G. F. Giudice, M. A. Luty, 
H. Murayama and R. Rattazzi, J. High Energy Phys. 12, 027 (1998). 

[21] G. Anderson et al, in New Directions for High Energy Physics, Snowmass 96, cd. by 
D. G. Cassel, L. Trindle Gcnnari and R. H. Sicmann; see also J. Amundson et al, in 
New Directions for High Energy Physics, Snowmass 96, ed. by D. G. Cassel, L. Trindle 
Gennari and R. H. Siemann. 

[22] Non-universal gaugino masses were first considered by, C. Hill, Phys. Lett. B135, 47 
(1984); Q. Shafi and C. Wetterich, Phys. Rev. Lett. 52, 875 (1984); J. Ellis, K. Enqvist, 

23 



D. Nanopoulos and K. Tamvakis, Phys. Lett. B155, 381 (1985); M. Drees, Phys. Lett. 
B158, 409 (1985). 

[23] C. S. Huang and Q.-S. Yan, Phys. Lett. B442, 209 (1998); C. Hamzaoui, M. Pospelov 
and M. Toharia, Phys. Rev. D 59, 095005 (1999); C. S. Huang, W. Liao and Q.-S. Yan 
Phys. Rev. D 59, 011701 (1999); S. Rai Choudhury and N. Gaur, Ref. [4]; K. Babu and 
C. Kolda, Ref. [5]; C. Bobeth et al. Ref. [8]; G. Isidori and A. Retico, Ref. [9]; Z. Xiong, 
J.M. Yang, Nucl. Phys. B628, 193 (2002); D. Demir, K. Ohve and M. Voloshin, hep- 
ph/0204119 (2002); A. J. Buras et al, Ref. [13]. 

[24] K. Babu and G. Kolda, hep-ph/0206310 (2002). 

[25] D. Pierce, J. Bagger, K. Matchev and R.-J. Zhang, Nucl Phys. B491, 3 (1997). 

[26] H. Baer, F. Paige, S. Protopopcscu and X. Tata, hep-ph/0001086 (2000). 

[27] D. Gronin-Hcnnessy et al. (Cleo Collaboration), Phys. Rev. Lett. 87, 251808 (2001); 

R. Barate et al. (Aleph GoUaboration) , Phys. Lett. B429, 169 (1998); K. Abe et al. 

(Belle GoUaboration), Phys. Lett. B511, 151 (2001). 
[28] H. N. Brown et al, (E821 GoUaboration), Phys. Rev. Lett. 86, 2227 (2001). 
[29] Joint HEP2 Supersymmetry Working Group, Combined Chargino Results, up to 208 

GeV, 

http://alephwww.cern.ch/lepsusy/www/ inos_moriond01/charginos .pub. html . ; 
Combined LEP Selectron/Smuon/Stau Results, 183-208 GeV, 

http : //alephwww . cern . ch/~ganis/SUSYWG/SLEP/sleptons_2k01 . html . 

[30] LEP Higgs Working Group GoUaboration, hep-cx/0107030. 

[31] H. Baer and M. Brhlik, Phys. Rev. D 55, 3201 (1997); H. Baer, M. Brhlik, D. Gastano 
and X. Tata, Phys. Rev. D 58, 015007 (1998). 

[32] S. Abel et al (SUGRA Working Group GoUaboration), hep-ph/0003154 (2000). 

[33] A. Gzarnecki and W. Marciano, Phys. Rev. D 64, 013014 (2001). This uses the analysis 
of M. Davier and A. Hooker, Phys. Lett. B435, 427 (1998), which in turn uses r decay 
data to reduce the error on the hadronic vacuum polarization. See also S. Narison, Phys. 
Lett. 513, 53 (2001). For a dicussion of the correction to the sign of the hadronic light 
by light contribution that reduces the significance of the originally reported result, see 
M. Knecht and A. Nyffeler, Phys. Rev. D 65, 073034 (2002); M. Knecht, A. Nyffeler, 
M. Perrottet and E. De Rafael, Phys. Rev. Lett. 88, 071802 (2002); M. Hayakawa and 
T. Kinoshita, hep-ph/0112102 (2001); J. Bijnens, E. Pallante, and J. Prades, Nucl. 
Phys. B626, 410 (2002); I. Blokland, A. Gzarnecki and K. Melnikov, Phys. Rev. Lett. 
88, 071803 (2002). 

[34] K. Melnikov, Int. J. Mod. Phys. A16, 4591 (2001). This uses the analysis of 
F. Jegerlehner, hep-ph/0104304. Melnikov's updated analysis of the SM value of Sa^ 
was presented at the High Energy Physics Seminar, University of Hawaii, March 2002. 

[35] E. Polycarpo, LHCb-2002-027 (unpublished); for an older analysis that also includes a 
discussion of the capabilities of ATLAS and GMS, see P. BaU, hep-ph/0003238 (2000). 
We note that conclusions from Ball's analysis are significantly more optimistic than 
those from Polycarpo. 

[36] H. Baer et al hep-ph/0008070 (2000). 

[37] H. Baer, M. Drees and X. Tata, Phys. Rev. D41, 3414 (1990); J. Ellis, G. Ridolfi and 

F. Zwirner, Phys. Lett. B237, 423 (1990). 
[38] H. Baer, P. Mercadante, X. Tata and Y. Wang, Phys. Rev. D 60, 055001 (1999). 



24 



[39] H. Baer, P. Mercadante, X. Tata and Y. Wang, Phys. Rev. D 62, 095007 (2000). 

[40] H. Baer, J. K. Mizukoshi and X. Tata, Phys. Lett. B488, 367 (2000). 

[41] H. Baer, M. Diaz, P. Quintana and X. Tata, J. High Energy Phys. 04, 016 (2000). 

[42] G. Anderson, H. Baer, C. H. Chen and X. Tata, Phys. Rev. D 61, 095005 (2000). 

[43] H. Baer, C. Balazs, J. Ferrandis and X. Tata, Phys. Rev. D 64, 035004 (2001). 

[44] D. E. Kaplan, G. D. Kribs and M. Schmaltz, Phys. Rev. D 62, 035010 (2000); Z. Chacko, 
M. A. Luty, A. E. Nelson and E. Ponton, J. High Energy Phys. 01, 003 (2000). Sec also, 
M. Schmaltz and W. Skiba, Phys. Rev. D 62, 095005 (2000) and Phys. Rev. D 62, 
095004 (2000). 

[45] H. Baer, A. Belyaev, T. Krupovnickas and X. Tata, Phys. Rev. D 65, 075024 (2002). 

[46] H. Baer et al, Phys. Rev. D 61, 111701 (2000); and D 63, 015007 (2000); H. Baer and 
J. Ferrandis, Phys. Rev. Lett. 87, 211803 (2001); T. Blazck, R. Dcrmisck and S. Raby, 
Phys. Rev. Lett. 88, 111804 (2002) and Phys. Rev. D 65, 115004 (2002), and references 
therein. 

[47] Y. Fukuda et al. Phys. Rev. Lett. 82, 2644 (1999) and 85, 3999 (2000). 



25 



TABLES 
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Mz 






M3 


M2 


Ml 


Ms 
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~ 6 


~ 6 
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200 
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10 


~ 6 


^ 4 


~ 10 



TABLE I. Relative gaugino masses at Mqut and Mz in the four possible irreducible represen- 
tations that F$ could transform as. 
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mSUGRA, = 300 GeV, m^^ = 300 GeV 
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FIG. 1. Branching fractions for the decays Bs l^^f^^ Siiid t^t~ in the mSUGRA 

model with mo = 300 GeV and mi/2 = 300 GeV, for > (sohd) and < (dashed). In frame 
a) we show the branching fractions versus tan (3 for = 0) while in frame h) we show these versus 
Aq for tan (3 = 40. Also shown is the branching fraction for the decay b ^ sj which is to be read 
off using the scale on the right. The squares mark the limits of the experimentally allowed regions 
as discussed in the text. 
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mSUGRA, = GeV 
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FIG. 2. Contours of constant branching fraction for the decays Bg 1^^ fJ-^ (sohd) and 
Bd T'^T~ (dashed) marked on the curves. The results are shown in the niQ — mi/2 plane of 
the mSUGRA model with Aq = 0, for a) tan/? = 50,/i > 0, and b) tan/? = 45,/i < 0. The dotted 
contour labelled 10~^ shows the result for B{Bs if the gluino loop contribution discussed 

in the text is set to zero. The dark-shaded region is excluded by theoretical contraints on EWSB. 
In the slant-hatched region, Zi is not the LSP. The region covered by open circles is excluded 
by experimental constraints on sparticle masses or on niA, as discussed in the text. Below the 
dot-dashed curve labelled h, ruh < 113 GeV. In frame a) below the solid curve labelled h 57, 
B{h — S7) < 2 X 10~^. Below the corresponding curve in frame 6), B{h — 57) > 5 x 10~^. 
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FIG. 3. Contours in the mp — rni/2 plane of the mSUGRA model with = where the 
branching fraction for Bg (B^ —>■ r+r") is 10~^ (10~^) for several values of tan/3. Shaded 

regions are as in the previous figure. In frame a) we take /x > 0, while in frame b) we take /x < 0. 
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mSUGRA, -3 < < 3 m 
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FIG. 4. The branching fraction for the decay Bg — > ^ scan of the mSUGRA parameter 

space over the range mentioned in the text. We show results versus mi/2 for a) tan/3 = 55,/x > 0, 
b) tan/3 = 45, /i < 0, and c) tan /3 = 40, /i < 0. Each cross (dot) denotes a model where B{h ^ sj) 
lies within (outside) the range (2 — 5) x 10~^. 
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FIG. 5. Contours of constant branching fraction for the decays Bg 1^^ 1^" (sohd) and 
Bd T~^T~ (dashed) with values marked on the curves in the A — tan (3 plane of the minimal 
GMSB model. We take the messenger scale M = 3A and ns = 1 in this figure and show results 
for a) /z > 0, and 6) /x < 0. In both frames, the dark-shaded region is excluded by theoretical 
considerations discussed in the text. In the slant-hatched region < 76 GeV, in the region 
covered by open circles mg^ < 100 GeV, and in the region covered by triangles, < 100 GeV. 
The open-circle region extends all the way to the smallest values of A (as does the slant-hatched 
region) but has been terminated at A = 30 TeV for clarity. To the right of the contour labelled 
6 — S7 in frame 6), B{h ^ 57) is in the range (2 — 5) x 10~^. It is in this range all over the plane 
in frame a). 
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GMSB, M = 3 A, tan P = 40 
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FIG. 6. The dependence of the branching fractions for the decays Bg /^^/^^ and t~^t~ 
on the parameter n^. We plot the branching fractions versus a sparticle mass (we choose rrig) 
rather than the theoretical parameter A for reasons discussed in the text. We show our results for 
a) /X > 0, and 6) /x < 0. Other parameters are as shown on the figure. 
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FIG. 7. Branching fractions for tlic decays Bs 1^-^ 1^-^ and t^t^ within the framework 

of the minimal anomaly-mediated SUSY breaking model. We show the branching fractions versus 
tan/3 for a) mo = 400 GeV and b) mo = 1 TeV, and other parameters as labelled on the figure. 
Also shown is the branching fraction for the decay b ^ sj which should be read using the scale on 
the right. 
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AMSB, tan |3 = 42, |i< 
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FIG. 8. Contours of branching fraction for the decays Bg —>■ (solid) and — r+r~ 

(dashed) in the mo — ^3/2 plane of the mAMSB model. The values of the branching fractions 
are mentioned in the text. The dark-shaded region is excluded by theoretical considerations. 
Along the squares and triangles running along the boundary of the upper dark-shaded region, nih 
and tha respectively fall below their experimental bound. In the region covered by open circles, 
< 100 GeV. Finally, in the region covered by dots close to where the squares and triangles 
are, B{b ^ 57) > 5 X 10-^, while over the unshaded region it is in the range (2 — 5) x 10 . 



34 



non-univ, m. = 400 GeV, A. = GeV, tan P = 40 
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FIG. 9. The branching fraction for the decay Bg /^^/^^ versus the GUT scale gluino mass 
parameter M3 in SU (5) supergravity models with non-universal gaugino masses. As discussed in 
the text, the models are characterized by the gauge transformation property of the auxiliary field 
F$ whose vev breaks SUSY and gives rise to gaugino masses. We show the branching fraction for 
the case where this field transforms as a singlet (dotted), a 24 dimensional representation (solid) 
and a 200 dimensional representation (dashed) of S\J (5) for a) > 0, and b) /j, < 0. 
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non-univ, tan (3 = 55, -3 hIq < < 3 m^, |i > 




-6 : 
10 s 



^^■ ■ ■ ■■■■■■■ ■iiiiii i ii i l i iii i 



I c)200DR 




^^^^^ ^^^^^ llllllgllll 
■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ 



200 300 400 500 600 700 800 

Ml (GeV) 

FIG. 10. The branching fraction for the decay Bs jJ'^ for a scan of the parameter space 
of SU{5) supergravity models with non-universal gaugino masses over the range mentioned in the 
text. We take /x > 0. We show results versus the GUT scale gluino mass parameter Mg for the 
a) 24 model, b) 75 model, and c) 200 model discussed in the text. Each cross (dot) denotes a 
model where B{b —>■ 57) lies within (outside) the range (2 — 5) x 10~^. 
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FIG. 11. The same as Fig. 10, except for fi <0. 
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